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Abstract 

The systematization of the purely Lagrangean approach to constrained 
systems in the form of an algorithm involves the iterative construction 
of a generalized Hessian matrix W taking a rectangular form. This 
Hessian will exhibit as many left zero-modes as there are Lagrangean 
constraints in the theory. We apply this approach to a general La- 
grangean in the first order formulation and show how the seemingly 
overdetermined set of equations is solved for the velocities by suitably 
extending W to a rectangular matrix. As a byproduct we thereby 
demonstrate the equivalence of the Lagrangean approach to the tra- 
ditional Dirac-approach. By making use of this equivalence we show 
that a recently proposed symplectic algorithm does not necessarily 
reproduce the full constraint structure of the traditional Dirac algo- 
rithm. 
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1 Introduction 



A number of algorithms have been developed over the past years for treat- 
ing constrained Hamiltonian systems. Perhaps the most familiar one to the 
physicist community is the one developed by Dirac Although very elegant 
and powerful in its algebraic structure, this algorithm has been critizized for 
being based on the existence of so called "primary constraints", which are 
a purely phase-space artefact, and have no counterpart on the Lagrangean 
level. Faddeev and Jackiw [|] have thus proposed an alternative method 
based on a first order Lagrangean (symplectic) formulation, avoiding the in- 
troduction of primariy constraints. Furthermore, the local symmetries of the 
Hamiltonian, as generated by the so called "first-class" constraints in Dirac's 
terminology, turn out to be larger than those of the Lagrangean This has 
led to a renewed interest in the problem of deducing the local symmetries 
of a Lagrangean from the Hamiltonian formalism j3], and in particular to 
a revival of the "Lagrangean approach", and the "symplectic approach" to 
constrained systems H, |5|, || [7|, ||, |[|. Of all three methods, the Lagrangean 
algorithm is actually the most pedestrian one, with a solid mathematical 
basis. The symplectic algorithm, on the other hand, as developed in a series 
of elegant papers in refs. j§], lacks a rigorous mathematical justification, and 
can lead, as we shall demonstrate, to an incomplete solution of the prob- 
lem. In order to establish the relation between the three formalisms, we shall 
thus take as our starting point the Lagrangean approach as applied to first 
order Lagrangeans, in order to allow for a comparison with the symplectic 
approach. As shown in section 2, the Lagrangean algorithm leads to a larger 
set of equations than the number of unknown velocities to be solved for. 
This is reflected in the fact that the generalized Hessian which implements 
the algorithm, is a rectangular matrix possessing as many left zero-modes, 
as there are Lagrangean constraints hidden in the Euler-Lagrange equations. 
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We show that these zero-modes are of such a form, that they permit the solu- 
tion of the equations of motion in terms of the inverse of a quadratic matrix, 
whose elements are just the Poisson brackets of the Hamiltonian constraints 
- including the primary constraints. We thereby establish the equivalence 
with Dirac's algorithm. In section 3 we then consider the simple example 
of the particle motion on a hypersphere and thereby demonstrate that the 
symplectic algorithm of ref. Q is not always equivalent to the Dirac and 
Lagrangean approach. We conclude this section by discussing the general 
condition under which this symplectic algorithm fails. Section 5 summarizes 
our findings. 



2 The Lagrangean algorithm 

Given a second order Lagrangean, one can always find an equivalent first 
order Lagrangean of the form 

L(Q,Q)=a a (Q)Q a -V(Q) (1) 

where Q stands for n degrees of freeedom Q a , a = 1, 2 • •-, n. The correspond- 
ing Euler-Lagrange equations read 



= (2) 

where the matrix is defined by 

W^(Q) = ^~ (3) 



* Wj dQ a dQ p 



Let r be the rank of the matrix W^°\ Then there exist n — r zero 
modes of which we denote by u^(a) , a — 1 • • • n — Tq. Multiplying 

equations (0) from the left with these zero modes, we are led to the zero-level 
Lagrangean constraints 

r)V 

^ 0) = E^ 0) H^ = 0, a = l,..-,n . (4) 



Some of these constraints may vanish identically. The remaining ones we 
denote by <p$. The corresponding zero modes u^°\ao) we refer to as "non- 
trivial" . 

In general there are further constraints hidden in equations (|J). In order 
to unravel them, we implement their conservation by adjoining their time 
derivatives 

to the equations (§). This leads to the following enlarged set of equations 

W Mp{Q)Qp = Kj&iQ) (6) 
where are now the elements of a rectangular matrix 



(7) 



with 



and 



U (0, _ dtp® 

7x77 {b) 




= (9) 



where 

= W(Q) (10) 

We now look for "non-trivial" zero modes («W(oi),Oi = 1,- ■ -.ni) of iyW, 
and repeat the steps above, adjoining the time derivative of any new con- 
straints to the equations of motion (H). Repeating this algorithm, the itera- 
tive process terminates after L steps, when no new constraints are generated. 

Denote the full set of constraints generated by the algorithm collectively 
by {ip a }, a = 1,---,N. Denote further the set {a, a} collectively by {A}. 
The final set of equations can then be written in the form 

W A pQp = K A (11) 



where 



W. 



A0 



\ M al3 



with 



and 



M. 



a/3 



dQ/s 



K(o) 



\ 







:i2) 



(13) 



(14) 



Denoting by u(a) the left zero-modes of the matrix Wa/3, the constraints are 
given by <p a = u(a) ■ K = 0. 

Equations ( |TTj) represent n + N equations for the n velocities {Q a }- In 
general such a set of equations would be overdetermined and admit no non- 
trivial solution. Since the additional iV equations were however generated by 
a self-consistent algorithm from the original Euler-Lagrange equations, the 



equations (|TTJ) do in fact admit a non-trivial solution. In the following we 
shall assume the first order Lagrangean ([]]) to describe a purely second class 
system in the Dirac terminology. In that case we have the following 



Assertion: 

The unique solution to fllTD for the velocities is given by 



A - p-i/r (0) 



(15) 



where F is the inverse of the matrix F obtained by extending the rectan- 



gular matrix W defined in (|12[ ) to the antisymmetric square matrix 



<AB 



-Ml 



ab 







(16) 



with M a(] defined in (0), and Mj b = M ba . 



Proof of Assertion: 
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Consider an enlarged space on which the square matrix (|16|) is to act (we 
streamline the notation in a self-evident way), 



£a := (Qa,Pa) 



(17) 



and the following equations: 



FabZb = K A 



As we shall prove further below, det F^O for a second class system. Hence 
we can solve these equations for the velocities 



U = F- A lK B 



(19) 



We write the inverse matrix F 1 in the form 



F ab 



( w<f>) 



aii 



■ML 



\ M af3 u ab 



Then F X F = 1 and FF 1 = 1 respectively imply, 



on rv ~fl3 ~ lv± ac^cf3 — > 

Wj^Ml = 
M ai W$ + u ac M cP = 
M ai M? b = -5 ab , 



Ja/3 



(20) 



(21) 
(22) 



and 



w<$w®-mLm, 



M m W$ = 



c/3 — u a/3 



W^Ml + MLu cb = 



Q"7 yb 



M n ^Ml, = -5 ab . 



(23) 
(24) 



Consider eqs. ( |i8|) which, written out explictely, read 



Ml b p b = tf(°> 



(25) 



MapQp = (26) 

From (P~5|) we see that the last equation states that ip a = 0, where ip a = 
are the constraints hidden in the equations of motion (0). Because of this, 
requiring their presistance in time implies that the second term on the LHS 
of (p5[) must vanish. Making use of (|22|) this in turn implies that p a = for 
all a. f\ Setting p a = in (|19"D , we have from (j20|). 



4* = W$tf< 0) (27) 

= M a pKf (28) 

Eq. ( f2~8|) is just the statement that = 0. To see this we notice that 
according to (0), the vectors 

u A (a) := (M aj ,u ac ) (29) 



are just the left zero modes of the matrix ([12]) . Hence 



M a/3 Kf = « A (a)i^ = • (30) 

As we now show, eqs. (|19"D are nothing but the Hamilton equations of motion 
derived from the so-called extended Hamiltonian. By making contact with 
the Hamiltonian formalism, we will prove that i) F is an invertible matrix, 
and ii) the solutions to ( p5|) imply that p a = 0, as was claimed above. This, 
at the same time, will prove the uniqueness of the solution. 

From the Dirac point of view, the symplectic Lagrangean (0) descibes a 
system with a primary constraint for every coordinate Q a : 

<Pa '■= P a - a a (Q) = , a = l,---,n (31) 

where P a are the momenta canonically conjugate to the coordinates Q a . 
Since the Lagrangean (HD is first order in the time derivatives, the corre- 
sponding canonical Hamiltonian H c is just given by the potential V, 

H c = V(Q) (32) 



3 An alternative proof that this must indeed be the case will be given further below, 
where we make contact with the Hamiltonian formalism. 



and hence does not depend on the momenta. The dependence on the mo- 
menta enters only in the total Hamiltonian via the primary constraints: 

H T {Q, P) = V(Q) + £ v a <j> a (Q, P) • (33) 

a 

The Dirac algorithm will in general lead to secondary constraints, which we 
label by a latin index: ip a = 0. It is easy to see that they are identical with 
the constraints generated by the Lagrangean algorithm. Thus consider the 
persistance equations for the primary constraints (j) a : 

h t } = {<p a , v} + Y,{^ = o . 



From (j3ll) and (§) we see that {4> a ,4>p} = W^){Q), so that the above eqia- 



tions read, 

W$jlv fi = K<® . (34) 

Multiplying this equation with the left-zero modes of we arrive at the 
level-zero Lagrangean constraints @, which are only functions of Q. Re- 

a/3 



quiring their persistance in time as generated by Ht yields Mjfjvp = 0, and 



adjoining these equations to ( 54]) 



<U = Kg • 

By taking appropriate linear combinations of these equations, new con- 
straints may be generated which are functions of only the QaS. This just 
corresponds to looking for left-zero modes of W^. The new constraints are 
thus identical with those derived in the Lagrangean approach at level "one" . 
Proceeding in this way it is easy to see that the secondary constraints gen- 
erated by the Dirac algorithm applied to Ht{Q, P) are identical with the 
constraints {ip a = 0} generated by the Lagrangean algorithm. 

We now go over to the extended Hamiltonian by including the secondary 
constraints with their respective Lagrange multipliers v a , 

H t ^H e = H c + J2 X b^b , (35) 

B 



where 

:= (<Pa, fa) , ^A-=(v a ,V a ). 

The Hamilton equations of motion for the coordinates Q a associated with 
the extended Hamiltonian He, read 

Q a = {Q a ,H E } = v a , (36) 

P a = {P a , H E } = -K® + Vp d a a p - v b d a cp b (37) 
Consistency with the persistance in time of the primary constraints requires 

0a = Pa — = (38) 

One readily verifies from ([37]) that this leads to 

W^v p - M T ah v b = JtJO . (39) 
On the other hand, persistance of the secondary constraints tp a leads to 

M bf ,Vf, = (40) 



Upon making use of (|3q) , we thus retrieved equations (P5|) and (|26|) if we 
identify v a with p a . Hence in the Hamiltonian formalism these equations are 
merely the persistance equations of the primary and secondary constraints, 
which can be compactly written in the Hamiltonian form 

{n A , h c } + Y,{^a, n B }\ B = o . (4i) 

B 



We now recognize that the matrix elements of F in ( fig) are given by 

Fab = {&a, &b} ■ 

Since the constraints have been assumed to be second class, this matrix is 
invertible. Noting further that 

{n A ,H c } = -K A , 



it follows from ( 4"I|) that 



A a = F A gK B ■ 



(42) 



With Q, 



a 



v a , these equations are nothing but fll9D , with p a identified with 



v, 



a- 



To prove the equivalence of (p9|) and (|40"|) with the original set of equations 
(p!l|), we must still show that eqs. ( |4*2| ) imply that v a = 0. To this effect we 
recall that the secondary constraints tp a = have actually been generated by 
the total Hamiltonian Ht from the persistance equations 



or equivalently M^ b Vb = 0, which, upon making use of fl22"|), implies v a = 0. 
This completes the proof of our assertion. 

Concluding this section we have therefore shown the full equivalence of 
the Lagrangean and Hamiltonian approach to the theory described by the 
first order Lagrangean ([!]). Any other approach must therefore reproduce 
the constrained structure of the Lagrangean approach. An alternative (sym- 
plectic) algorithm for unravelling the constrained structure was proposed in 
ref. [4]. In the following section we will show that the symplectic algorithm 
does not necessarily generate the correct constrained structure. 

3 The symplectic algorithm 

In the following we first illustrate in terms of a simple example, an alternative 
algorithm for generating the constraints, as proposed in ref. [|J. We refer to 
it as the "symplectic algorithm" . 



Consistency with (f4i"l) therefore requires that 
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3.1 Particle on a Hypersphere 

The following (second order) Lagrangean is referred to as describing the 
linear sigma model in Quantum Mechanics: 

where q — (qi, ■ • • , q n ). The equivalent symplectic Lagrangean reads, 



with 



y(0) 



-Kf - 1) + \f 



The Lagrangean is of the form 



LW=a a (Q)Q a -V(°\Q) 



with 

Qa ■= (q,p t A) , a« := (jP, 0, 0) 
The equations of motion are of the form with 



W {0) 



0-10 
10 

T T 



and 



V 



The matrix has one "zero-level" zero mode: 



«<?> := (0,0,1) 



implying the constraint 



<p<® = -u^K^ = g 2 - 1 = 
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This constraint will necessarily coincide with that of the Lagrangean ap- 
proach at the zero'th level. 

In the symplectic algorithm the time derivative of the constraint fl5"0|) is 
however added in the (partially integrated) form —p^cp^ to the zero- level 
Lagrangean (46), [] to yield the first level Lagrangean 

L« = L<°> - p(°yo) , (51) 

where p^ is a new dynamical variable. Correspondingly we define the ex- 
tended set of coordinates 

can be written in the form 

LV = ami-VV>(Q), 

where 

ag:= (p,0,0,-^°)). 
The corresponding Euler-Lagrange equations read 

where the "first-level" symplectic square matrix is given by, 



and 



( 


-1 




1 










T 




Q T 




/ 


-2Ag 






P 






(<?-! 




V 






/ 



(52) 



(53) 



4 In ref. Q the term proportional to the Lagrange multiplier A in L' ' has been absorbed 
into the term proportional to p. 
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has two (level-one) zero modes, 

^(1):= (0,0,1,0) 

u£(2):=(0,2g,0,-l) 

The first zero mode reproduces the constraint (f — 1 =0. The second zero 
mode yields the new constraint 

V9 (1) = u {l \2)-K {l) = 2p-q 

As one easily verifies, these constraints are identical with those obtained in 
the Lagrangean algorithm described in section 1, at this level. 

According to the symplectic algorithm we now define the second level 
Lagrangean by adding the new constraint in the form 



L (2) = L (0)_ -(0)^(0) (1)^(1 



or 



with 



and 



r(2) _ „( 2 )t( 2 ) 



„(2) ._ 
a A 2 ■- 



(p : 0A-(f-l),-2p-q) 



For the corresponding symplectic matrix one obtains 



i?(2) 






-1 


— * 




-2q 


-2p 


1 





— * 




— * 




-2q 


T 


0^ 











2f 


T 











2f 


2f 












(54) 



/ 

As one readily checks, this matrix has only one zero mode u^] '■= (0 T , T , 1,0,0) 
which, however, just reproduces the constraint ip^ = 0. Hence the sym- 
plectic algorithm terminates at this point, leaving us with a non-invertible 
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matrix. Q On the other hand, one readily checks that the standard La- 
grangean (or equivalently, Dirac) algorithm generates not only the constraints 
<f — 1 = 0, p- q = 0, but also one futher constraint 2\q* + fp = 0. Indeed, in 



the Lagragian algorithm, in (5l) is replaced by the rectangular matrix 






-1 


1 





6 T 


T 



2<f 1 



y 2f 2? 
which is seen to possess the three zero modes, 



(55) 




1 




u^(2) 



f o N 

2q 


-1 

V J 



< 2t \ 



-2p 


1 



V 



(56) 



which in addition to the constraints (p^ = 0, ip^ 1 ' = 0, imply a new constraint 

if^ := 2\q> + f = . (57) 



Hence we are taken to a third level with the corresponding enlarged matrix 
given by 



/ 



W {3) 





1 

T 
2<f 
2^ 





5* 



2q 



7fT 



\ 



4A«f 2f 2? 



(5* 



5 In ref. j| the Lagrange multiplier A was absorbed into the dynamical variable p(°\ 
Theirby the information about A was lost, and the resulting matrix at level 2 became 
invertible. 
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As one readily checks, has no new zero modes. Hence the algorithm 

terminates at this point. Notice that the extension of this matrix to a square 
matrix as discussed in section 2 results in an invertible matrix, reflecting a 
second class system. 

We see that the symplectic algorithm fails to generate the correct set of 
constraints known to be present for the model in question. In fact, from the 
point of view of the second order Lagrangean formulation there exists just 
one primary constraint <fi = p\ = 0, where p\ is the momentum conjugate 
to the variable A, and the total Hamiltonian correspondingly reads, Ht = 
\fr 2 — X((f — 1) + vp x . As one readily checks, the last constraint (|57|) just 
serves to fix the Lagrange multiplier v in Ht to v = 0. Only at this final 
stage the second class nature of the model in question becomes evident. If 
we stop at level two, v remains arbitrary, as expressed by the zero column in 
(H) and (H). 

3.2 When does the symplectic algorithm fail? 

We now examine in general terms at which point the symplectic algorithm 
begins to fail. To this end we examine what the symplectic algorithm de- 
scribed above corresponds to on Hamiltonian level. Let l/ -* be of the form 
([]]), with ( |3TD the corresponding primary constraints. At the I + l'th level, 
the symplectic algorithm leads to a Lagrangean of the form (we streamline 
the notation) 

L (m ) = L (o)_ E/Wg)i 

where (p av a e = 1, • ■ • , ng denote all the constraints generated by the iterative 
procedure up to level I. The corresponding total Hamiltonian reads, 

Ht +X = #r 0) + Ki^a t ■ 
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Here {<f) ae } denote the corresponding set of primary constraints associated 
with {p a J, 

<fia e = + ( Pa e (Q) . 

where P ae are the momenta conjugate to the dynamical variables p ae and 

with (p a the primary constraints (0), associated with the original Lagrangean 
L^°\ Hence in the symplectic algorithm described above, the total Hamilto- 
nian is modified at each level. Clearly the Euler-Lagrange equations derived 
from L( £+1 ) and the Hamilton equations of motion following from H^ +1 > de- 
scribe the same dynamics. 

Conservation of all the primary constraints requires, 

{<f> at ,H$ +1) } = £{^,^w + £{<^,&ja 6/ = o. ( 59 ) 

Let $A e stand for 

with Xa,, the corresponding set of Lagrange multipliers: 

^A e '■= (v a , X ae ) 

Then we may write (fj3) in the compact form 

£{Q^,0 B JA^ = i^, (60) 

Bi 



where 



K = (W, 0) 



with an Ng, = n + ng - component null- vector. 

One readily checks that {f^,0# f } is identical with FA e B e in ( P>D at the 
£'th level. Furthermore, with the identification of v a and A a< , with Q Q and 
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p ae via the Hamilton equations of motion, 

Qa = {Qai Hrp j- = V a , 
Pag — \Pa n tl T ) — A a( 

we see that the persistance equations (|59D are just the equations of motion 
obtained from L( £+1 ) in the symplectic approach. 

Within the Hamiltonian formalism, the search for zero modes of F at 
level £ now corresponds to seeking linear combinations of all the primaries, 
ua^Ao sucn th & t 

= 0. (61) 

From (p0| ) we see that these equations imply linearly independent (non- 
trivial) constraints, which we denote by 

Pan = u Ae (ai)K Ae = . 

Of the conditions flBTD , only those with Be = /3, 

Y, u AMd{ ( $>A l Ap} = ® (62) 

At 

are contained in the Lagrangean (and hence traditional Dirac) approach. Let 
u(ag) be solutions of (|62|) . From (pi]), with B e = hi we see that the symplectic 
approach thus implies the additional restrictions 

u a (a e )^ = 0, (63) 

for the zero modes, which are not contained in the Lagrangean algorithm. 
Hence we have a mismatch between the symplectic and Lagrangean algo- 
rithm, once the latter condition is not satisfied in the iterative process, and 
the constraint structure becomes inequivalent for the two algorithms. This 
is the main point of this paper. 

6 This subtle point has been missed in ref. [0 
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Let us exemplify this for the case of the "particle on a hypersphere" . At 
the second level the Lagrangean algorithm leads to the zero modes (0). We 
verify that at the zero'th and first level of the iterative process the condition 
(|33|) is still verified, whereas this is not the case for the second level zero 
mode «( 2 )(3) in (|56|), since 

This explains why the algorithm stops before generating one further con- 
straint, ip& = 0, eq. 

It is instructive to further ellucidate the meaning of this finding. Going 
through the iterative procedure on Hamiltonian level (found above to be 
equivalent to the symplectic algorithm), we arrive after the second iterative 
step at the Hamiltonian 

n 

H® = v^(q, P ,x)+T t K( p qi -Pi)+^ i Pp i ) 

1=1 

+ vp x + A x (Pi +f-l) + A 2 (P 2 + 2q ■ p) 

Conservation in time of the primaries now merely serves to fix all the La- 
grange multipliers A, Ai, A 2 , A 2 , and leads to: 

v p + 2X x q + 2X 2 p = 2Xq 
v q + 2X 2 q = p 
Vq-p + v p - q = 
Vq-q = , 

as well as the constraint (f — 1 =0. These equations may be solved for Ai 
and A 2 , 

Ai = l -{f? + 2X-2p-q) 

A 2 = ^{p-q), 

and hence vor v qi and v Pi , leaving v undetermined. Hence <f — 1 = is the 
only constraint (as reflected by the zero column in ([54]) and ([55])), unless we 
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set Ai = A2 = 0. This just corresponds to working with the total Hamiltonian 
HP. In that case the algorithm does not terminate, but rather generates 
one further constraint, p 2 + 2 A = 0, whose time independence will finally fix 
also v to vanish. 

4 Conclusion 

In this paper we have examined the interrelation between three different al- 
gorithms currently in use for unravelling the constrained structure of first 
order Langrangians. We have referred to these as the "Lagrangean" , "Dirac" 
and "symplectic" algorithms. Of these the first two rest on a solid founda- 
tion, and, as we have seen, there exists a one-to-one correspondance between 
these formalisms. In particular we have shown how to invert the seemingly 
overdetermined system of equations of the Lagrangean algorithm. As for the 
symplectic algorithm presented in refs. |4j], it does not always reproduce the 
correct set of constraints, as we have seen. In fact, we have shown for a 
general first order Lagrangean, under what conditions the algorithm fails to 
reproduce all of the constraints correctly. A concrete example has exemplified 
this. 
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